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ABSTRACT 

We place constraints on the redshift-averaged, effective value of the equation 
of state of dark energy, w, using only the absolute ages of Galactic stars and 
the observed position of the first peak in the angular power spectrum of the 
CMB. We find w < —0.8 at the 68% confidence level. If we further consider 
that w > — 1, this finding suggests that within our uncertainties, dark energy is 
indistinguishable from a classical vacuum energy term. 

We detect a correlation between the ages of the oldest galaxies and their 
redshift. This opens up the possibility of measuring w(z) by computing the 
relative ages of the oldest galaxies in the universe as a function of redshift, dz/dt. 
We show that this is a realistic possibility by computing dz/dt at z ~ from 
SDSS galaxies and obtain an independent estimate for the Hubble constant, 
H = 69 ± 12 km s _1 Mpc -1 . The small number of galaxies considered at z > 0.2 
does not yield, currently, a precise determination of w(z), but shows that the age- 
redshift relation is consistent with a Standard ACDM universe with w = —1. 
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1. Introduction 

There is now solid observational evidence that the universe is accelerating. The absolute 
ages of high-redshift galaxies (Dunlop et al. 1996; Spinrad et al. 1997) and the value of the 
Hubble constant H (Freedman et al. 2001), the Hubble diagram of Type la supernovae (Riess 
et al. 1998; Perlmutter et al. 1999), the angular power spectrum of the cosmic microwave 
background (CMB; Benoit et al. 2002; Ruhl et al. 2002; Goldstein et al. 2002) in combination 
with galaxy surveys (Efstathiou et al. 2002; Verde et al. 2002) and lensing (Wang et al. 2002), 
all indicate that the universe is flat and dominated at present by some form of dark energy 
with negative pressure. Strong confirmation for this comes from the newly released results 
from the Wilkinson Microwave Anisotropy Probe (WMAP) first year data (Spergel et al. 
2003). 

The equation of state of the dark energy, p = wp, expresses the ratio between the 
pressure, p, and the mass density, p, of the dark energy in terms of the parameter w (in units 
of c = 1). The value of w could either be constant, as in the case of a cosmo logical constant 
(w = — 1), or be time- dependent, as in the case of a rolling scalar field or "quintessence" 
(Peebles & Ratra 1988; Caldwell et al. 1998). Any such behavior would have far-reaching 
implications for particle physics. 

In this paper we use the ages of the oldest stellar populations between z = and z — 1.5 
to derive new and independent constraints on two key cosmological parameters, w and Hq. 
Both parameters are relatively poorly constrained by CMB data alone. For example, in order 
to get accurate constraints on w, Spergel et al. (2003) combined the WMAP data with six 
external data sets: two additional CMB data sets (CBI and ACBAR; Pearson et al. 2002; 
Kuo et al. 2002), large scale structure measurements from the 2dF galaxy redshift survey 
(2dFGRS; Percival et al. 2001; Verde et al. 2002), measurements of the Lyman a forest 
power spectrum (Croft et al. 2002; Gnedin & Hamilton 2002), the luminosity distance- 
redshift relation from Type la supernovae measurements (Garnavich et al. 1998; Riess et al. 
2001), and H from the Hubble Space Telescope (HST) Key Project (Freedman et al. 2001). 
Therefore it is particularly important to obtain accurate and independent measurements of 
cosmological parameters relying on different physics. 

In Section 2 we focus on w, summarizing the basis of our method and deriving tight 
constraints on the equation of state using the ages of globular clusters in combination with 
the location of the first acoustic peak in the CMB power spectrum. Then, in Section 3, we 
determine stellar ages for a sample of early-type galaxies between z = and z — 1.5 and use 
the age-redshift relationship defined by the old envelope to determine H . The prospects of 
using this method to derive a measurement of w(z) are also discussed. Section 4 summarizes 
and discusses the results. 
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Throughout the paper we adopt the following notation. The Hubble constant is H = 
100 h km s" 1 Mpc -1 and c is the speed of light. The density of baryons, matter, radiation, 
photons, neutrinos, and dark energy, at present, in critical units, are indicated as fib, fl m , 
^rad> ^7> ^v, an d respectively (fl ra d = fi 7 + fi„). The physical quantities, independent of 
h, are indicated with u; e.g., uo m = Vt m h 2 . When needed, a standard, flat, reference ACDM 
model with fl m = 0.27, h = 0.71, Qt>h 2 = 0.023, and w = — 1 is assumed (Spergel et al. 
2003). 

2. The value of w from stellar ages 

In this section we use the ages of the oldest globular clusters in combination with the 
location of the first acoustic peak of the CMB power spectrum (Page et al. 2003) to constrain 
the value of w. The method is briefly described in Section 2.1 and the results are given in 
Section 2.2. 



2.1. Method 

If the universe is assumed to be flat, the position of the first acoustic peak (£\ in the 
standard spherical harmonics notation) depends primarily on the age of the universe and 
on the effective value 7 of w (Caldwell et al. 1998; Hu et al. 2001; Knox et al. 2001). As 
noted by these authors, for a fixed w value, a change in the physical density parameter 
uj m = fl m h 2 that keeps the characteristic angular scale of the first acoustic peak fixed will 
also leave the age approximately unchanged. Thus an independent estimate of the absolute 
age of the universe at z = combined with a measurement of i\ yields an estimate of w 
largely independent of other cosmological parameters (e.g., Q m and h). 

This can be better understood by considering that for a constant w(z), the age of a flat 
universe is given by 

t = h q 1 f°(i + z)- 1 [n rad (i + Z ) A + vi m {\ + zf + n A (i + z) 3 ^] ~ 1/2 dz . (i) 

Jo 

The position of the first acoustic peak is fixed by the quantity 

9 A = r s (a dec )/d A (a dec ) , (2) 



7 I.e., the average value over redshift. We use w(z) to indicate when w is allowed to vary over time. 
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where ad ec is the scale factor at decoupling, r s (ad ec ) is the sound horizon at decoupling, and 
dA(ddec) is the angular diameter distance at decoupling. For a flat universe, 
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See Verde et al. (2003) for more details. We use equations 1 to 4 and fix 9 A to be consistent 
with the WMAP determination {t\ ~ 220; Page et al. 2003). Figure 1 shows the allowed 
region in the age-w plane obtained for Vl^ti 1 = 0.02, 0.05 < Q m < 0.4, and 0.5 < h < 0.9. It 
is clear from the plot that an independent and accurate age determination of the universe 
will provide a measurement of w. Thus by using only the WMAP observation of i\ and 
an independent estimate of the age of the universe, we can place constraints on w, largely 
independent of other cosmological parameters. 



2.2. Results from the ages of globular clusters 

The ages of the oldest globular clusters (GCs) provide a lower limit to the total age of 
the universe. Numerous star-forming galaxies have now been observed up to redshift z = 6.6 
(e.g., Kodaira et al. 2003). Since the oldest GCs contain the oldest stellar populations in 
galaxies, it is reasonable to assume that GCs too have formed at redshift z > 6. Furthermore, 
since small halos are easily destroyed after reionization (e.g., Bullock et al. 2001), GCs likely 
formed prior to the reionization epoch. Under this assumption, the WMAP recent estimate 
of the reionization redshift z rcion ~ 17 (Kogut et al. 2003) implies that the ages of the oldest 
globular cluster do not significantly underestimate the age of the universe. 

Globular cluster ages can reliably be estimated in several different ways: from the 
luminosity of the main-sequence turn-off (if the distance is known), from detached binaries, 
from the termination of the white-dwarf sequence, and from the luminosity function method 
which provides a joint constraint on the age and distance of a GC. Renzini et al. (1996) 
determine the distance to NGC 6752 by means of accurate photometry of GC white dwarfs 
and obtained an age of 14.5 ±1.5 Gyr. Using the same method they find 47 Tucanae to have 
an age of 13 ± 2.5 Gyr (Zoccali et al. 2001). Gratton et al. (1997), using the new Hipparcos 
distance scale to determine main-sequence turn-offs, estimate an age of II.8+I5 Gyr (95% 
confidence level) for the oldest Galactic GCs. Jimenez & Padoan (1998) apply the luminosity 
function method to M55 and obtain an age of 12.5 ±1.0 Gyr (95% confidence level). Hansen 
et al. (2002), using very deep HST images, find the white dwarf cooling sequence for M4 
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and estimate its age to be 12.7 ± 0.7 Gyr at the 95% confidence level. Using the theoretical 
age-mass relation and applying it to the detached eclipsing binary OGLEGC-17, Chaboyer 
& Krauss (2002) determine Omega Centauri to have an age of 11. 2 ±1.1 Gyr (95% confidence 
level). These age estimates are all remarkably consistent. Krauss & Chaboyer (2003) perform 
the most careful analysis to date of the effects of systematics in GC age determinations. They 
estimate the age of the oldest Galactic GCs using the main-sequence turn-off luminosity and 
evaluate the errors with Monte-Carlo techniques, paying careful attention to uncertainties 
in the distance, to systematics, and to model uncertainties. They find a best-fit age for the 
oldest Galactic GCs of 12.5tl'^ Gyr (95% confidence limits). We adopt their probability 
distribution for the oldest GC age, which we find can be accurately described by 



P(t) = ^T) eXP 



In [(t - T)/mf 
2o~ 2 



(5) 



where t denotes the age of the oldest GCs, A = 1.466, a = 0.25, T = 6.5 and m = 5.9. Since 
reionization took place at z ~ 17 (Kogut et al. 2003) and we assume that oldest GCs must 
have formed before then, for all reasonable cosmologies we only need add 0.2 — 0.3 Gyr to 
the GC ages to obtain an estimate of the age of the universe. We conservatively add 0.3 
Gyr 8 . 

In order to constrain w, we assume a flat universe and Monte Carlo simulate the distribu- 
tion of l\ subject to only weak constraints on the other cosmological parameters. For different 
values of w, we generate 10 5 models, randomly sampling the cosmological parameters Q m , 
Vt b h 2 and h with uniform priors, 0.5 < h < 0.9, 0.05 < Vt m < 0.45, and 0.01 < Vt b h 2 < 0.03. 
We then impose an age of the universe constraint by randomly sampling these models so 
that the age of the universe has a probability distribution whose shape is given by equation 
5, but offset by 0.3 Gyr. We use the publicly available code CMBFAST (Seljak & Zaldarriaga 
1996) to compute the resulting distribution of t\. This is shown in Figure 2 where the two 
solid lines are the 68% and 90% joint confidence levels. 

As expected, the age alone does not constrain w; there is a degeneracy between w and 
£ ± . If we now use the observed position of the first acoustic peak as recently measured, in 
a model independent way, from WMAP (Page et al. 2003), we can constrain w with high 
accuracy. The filled contours in Figure 2 are marginalized values for w at the 68% and 90% 
confidence levels. Thus we find w < —0.8 (w < —0.67) at the 68% (90%) confidence level. 
If we drop the assumption that the oldest GC must have formed before reionization, but 
still impose that they must have formed by z = 6, this constraint is slightly weakened. In 



8 The uncertainty in this correction is negligible compared to the uncertainty in the determination of GC 
ages. 
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this case the age distribution of equation 5 needs to be offset by < f Gyr and we obtain 
w < —0.5 at the 90% confidence level. We have also repeated the calculation using Gaussian 
priors for the cosmological parameters, h = 0.7 ± 0.2, Vt m = 0.25 ± 0.2, Vt b h 2 = 0.02 ± 0.01, 
and find similar, but slightly more stringent constraints on w. 

We stress here that this determination depends solely on the GC determination of the 
age of the universe and on the observed position of the first acoustic peak in the CMB 
power spectrum. Our constraint is slightly less stringent than that obtained by Spergel et 
al. (2003) from a joint likelihood analysis of WMAP with six external data sets ( WMAP 
+ CB1 + ACBAR + 2dFGRS + Lymana forest power spectrum + Type la supernovae + 
Ho constraint from the HST key project), but is tighter than the CMB-only {WMAP + 
CB1 + ACBAR) determination and comparable to the WMAP + ACBAR + CBI + HST 
constraint. 



3. Hq and w(z) from the ages of unresolved stellar populations at z — — 1.5 

3.1. Method 



We now turn our attention to the age of the oldest objects in the universe at z > as a 
function of redshift. As recently suggested by Jimenez & Loeb (2002), these can be used to 
determine H(z) = — Y+z^f an< ^ therefore to measure w(z) and H . The Hubble parameter 
is given by H 2 (z) = HQ[p tot (z) / p tot (0)], where p tot is the total density of the universe. For a 
flat universe composed of matter and dark energy, H(z) is related to the equation of state 
by 
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_ x dz 
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(6) 



where we have used the conservation equation for dark energy density, p\ = — 3H(z)[l + 
w{z)]p\. Thus a measurement of dz/dt is also a measurement of w(z); for further details, 
see Jimenez & Loeb (2002). 

To obtain dz/dt, Jimenez & Loeb (2002) proposed to use the old envelope of the age- 
redshift relation of E/S0 galaxies. We assume that E/S0 at different redshifts are drawn from 
the same parent population with the bulk of their stellar populations formed at relatively 
high redshift (e.g., Bower et al. 1992; Stanford et al. 1998). Then, at relatively low redshift, 
they are evolving passively and may be used as "cosmic chronometers" . 

In the rest of this section we will apply this method to a sample of spectra (described 
in §3.2) to identify the age-redshift relation defined by the old envelope (§3.3) and derive 
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an independent measurement of H (§3.4). We show that at the level of accuracy of our 
measurement, the use of the age-redshift relation allows us to apply this method even if a 
fraction of the stars in E/SO galaxies is assembled at z < 1 (e.g., Treu et al. 2002). We also 
briefly discuss the prospects of using this method to derive w(z). 

3.2. Sample and measurement 

In order to apply this method we need high-quality spectra of a sample of old stellar 
populations covering the largest possible range in redshift. To this aim we combine the 
following data sets: (i) the luminous red galaxy (LRG) sample from the Sloan Digital Sky 
Survey (SDSS) early data release (Eisenstein et al. 2001); (ii) the sample of field early-type 
galaxies from Treu et al. (1999, 2001, 2002, hereafter the Treu et al. sample); (iii) a sample of 
red galaxies in the galaxy cluster MS1054— 0321 at z — 0.83; and (iv) the two radio galaxies 
53W091 and 53W069 (Dunlop et al. 1996; Spinrad et al. 1997; Nolan et al. 2003, Dey et al., 
in preparation). 

The LRG sample is the most heterogenous, since it includes red galaxies in general, 
regardless of morphology. The Treu et al. sample contains high quality spectra with average 
S/N = 16 per A. The spectra in the Treu et al. sample have higher S/N than the typical 
LRG spectrum, but the wavelength coverage is better for the LRG sample (see example 
in Figure 3). To select potentially passively evolving galaxies, we include only galaxies 
with S/N > 9 per A and for which passively evolving synthetic stellar population models 
(Jimenez et al. 1998) provide a good fit to the continuum. This selection criterion discards 
dusty or star-forming objects. The spectra of seven red galaxies in the cluster MS1054— 0321 
were obtained by one of us (DS) on UT 11 March 2002 using the dual-beam Low Resolution 
Imaging Spectrometer (LRIS; Oke et al. 1995) at Keck Observatory. On the blue side, we 
used the 300 lines mm -1 grism (Abi aze = 5000 A; spectral resolution AAfwhm ~ 13.5 A) and 
on the red side, we used the the 400 lines mm -1 grating (A b i aze = 8500 A; spectral resolution 
AAfwhm ~ 8 A). These observations, which were obtained at the parallactic angle, totaled 
one hour of integration. In order to ensure no gaps in the spectra, half the data were obtained 
with a D560 dichroic and the remainder were obtained with a D680 dichroic. These spectra 
have a S/N of 8. We also add the red radio galaxies 53W091 and 53W069 for which accurate 
ages have already been determined (Dunlop et al. 1996; Spinrad et al. 1997; Nolan et al. 2003, 
Dey et al., in preparation). Examples of the spectra used are given in Figure 3. 

The age of the dominant stellar population in the galaxies is obtained by fitting single 
stellar population models (Jimenez et al. 1998, 2003a) to the observed spectrum. These 
models have been extensively tested (Kotoneva et al. 2002; Jimenez et al. 2003b) and have 



- 8- 



two free parameters: age and metallicity. By construction, star formation is assumed to 
occur in a single burst of duration much shorter (< 10%) than the current age of the galaxy 
(Jimenez et al. 1999). Therefore, the derived ages are single-burst equivalent ages. The 
best fit model is found by standard x 2 minimization. Then, stellar ages and their errors are 
obtained by marginalizing the resulting likelihood surface over metallicity. Some examples 
for the parameter confidence regions, at the 1 and 2a level, are shown in Figure 3. Note 
that for galaxies with high S/N, like the Treu et al. sample, the constraints in both age and 
metallicity are very tight and thus the so-called age-metallicity degeneracy is dramatically 
reduced and does not affect the age determination at a significant level. 

Before we can use the estimated ages to constrain H we have to estimate the effects of 
the potential systematic uncertainties related to adopting a single-burst stellar population 
model and a single-metallicity model. As far as the first uncertainty is concerned, we know 
that even massive field E/SO galaxies can experience some level of star formation activity 
at z > 0.2 (Menanteau et al. 2001; Treu et al. 2002; Willis et al. 2002). The presence of a 
young and luminous stellar population superimposed on a dominant old stellar population 
can bias the single-burst equivalent ages of the spectrum to be much younger than the mass- 
weighted average age. Although this is a serious concern for the interpretation of individual 
cases, our method relies on the old envelope of the age distribution at any given redshift. 
Therefore, the effects of this bias on the determination of H and w are strongly reduced, 
as we discuss further in § 3.2. Concerning single metallicity, the uncertainty connected to 
metallicity gradients is negligible with respect to other sources of uncertainty. In fact, our 
spectra typically cover 6 kpc (for the LRG) to 10 kpc (for the higher redshift samples). 
Given that typical metallicity gradients in E/SOs follow <ilog[Fe/H]/cflogr = —0.2 (Davies 
et al. 1993), metallicity only changes by 0.2 dex within this region. This is fairly small and 
the single-metallicity approximation affects our results well within the la value of the error 
in the recovered metallicity. 



3.3. An age-redshift relation and w(z) 

Figure 4 shows the derived single-burst equivalent ages of the galaxies as a function of 
redshift. The circles correspond to galaxies in the SDSS LRG sample, triangles to the Treu 
et al. sample, diamonds are galaxies in MS 1054— 0321, and crosses are 53W091 and 53W069. 
For clarity, the SDSS LRG points have not been plotted for ages < 7 Gyr. Typical errors 
on the ages of LRG galaxies are 10% and are not plotted. An age-redshift relation (i.e., an 
"edge" or "envelope" of the galaxy distribution in the age-redshift plane) is apparent from 
z = to z = 1.5. 
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Although an envelope can be identified, it is important to notice that there is a large 
spread of ages at any given redshift. In particular, the measured age for many galaxies makes 
them significantly younger than the age-redshift envelope. This is in part due to a real spread 
in the star formation histories of early-type galaxies, especially in the field, but mostly a 
result of using single-burst equivalent ages. Single-burst equivalent ages effectively measure 
a luminosity-weighted age. Therefore, even a relatively minor episode of star formation is 
sufficient to significantly underestimate the age of the bulk of (old) stars. Indeed, many of 
the galaxies in the Treu et al. sample show [O II] A3727 emission, consistent with minor 
ongoing star formation activity (Treu et al. 2002). Is this contamination enough to bias our 
results, or is the use of the relative ages of the old-envelope sufficient to reduce this bias to 
useful levels? 

To answer this question we consider a scenario where the bulk of stars in E/SO galaxies 
formed at high redshift, while a relatively small fraction (up to 10% in mass) of stars is 
formed at later times (z < 1.5). This scenario is approximately consistent with the observed 
evolution of the number density (Im et al. 2002), colors (Jimenez et al. 1999; Menanteau et al. 
2001), and spectral properties of early-type galaxies (Trager et al. 2000; Treu et al. 2002). 
Within this scenario, we use Monte Carlo simulations to estimate the effects of this bias on 
the age-redshift relation and ultimately on H and w(z) determinations. We generate galaxy 
spectra assuming that after the initial burst of star formation at z > 3 an extra 10% of the 
stellar mass is formed at a time t\ between redshift < z < 1.5 with uniform probability for 
t\. For simplicity, we use discrete values for the ages with spacing 0.1 Gyr and assume a flat 
VL m = 0.27, w — — 1 and h = 0.71 cosmology to convert time to redshift. The metallicity for 
a given galaxy is chosen to be uniformly distributed in the range 0.1 < Z/Z & < 4. Finally, 
we add random noise to the spectra so that S/N = 15 and determine the age with the same 
procedure used for the data. 

The results of our Monte Carlo simulations are shown in Figure 5. The top panel shows 
the simulated age-redshift distribution if only six early- types are found per redshift interval 9 . 
The edge is blurred and smeared due to recent star formation contamination and because 
of the error in the age determination; nevertheless, it is clearly recognizable. The galaxies 
observed within a few 100 Myr of bursts have very small single-burst equivalent ages, as 
observed for some galaxies in the LRG and Treu et al. sample, but they do not interfere with 
the location of the old envelope. 



9 We have drawn galaxies ages randomly from a grid that spans the age range 0.1 to 14 Gyr, with a 
spacing 0.1 Gyr. For our standard, baseline cosmology this age spacing translates to the redshift interval 
seen in Figure 4. 
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The lower panel shows the distribution for the case where 60 early-types are observed per 
redshift bin (this is less than the number of galaxies in the SDSS LRG sample). The solid line 
shows the theoretical position of the edge, which is remarkably close to observed edge. Note 
also that increasing the number of galaxies has improved the agreement between the observed 
and theoretical edge. Thus we conclude that small, recent bursts of star formation do not 
introduce a significant bias on the edge position; our primary uncertainty in determining 
w(z) to high redshift is the sparse size of the sample of distant, early-types considered here. 
We can now attempt to recover dz/dt from this simulation by looking at the shift in the 
edge position as a function of z (see § 3.4 and Jimenez & Loeb 2002) and thus deduce w(z). 
We recover w(z) — — 1 ± 0.2 for the whole redshift range. This shows that, with a sample 
similar to the LRG sample extended to higher redshift, we should be able to measure w(z) 
within a few 10%. This is similar to what was already concluded in Jimenez & Loeb (2002). 

Having established that the age-redshift relation for the old envelope is not significantly 
affected by the bias caused by later episodes of star formation, we now consider the obser- 
vational limits that we can impose on w(z). As a consistency check we note that the age 
at z = obtained here is in good agreement with GC ages (see §2) and that the Treu et 
al. galaxy ages agree with those from of SDSS LRG sample where they overlap. The solid 
line in Figure 4 corresponds to the age-redshift relation for our flat, reference ACDM model: 
Vl m = 0.27, h = 71, and w = —1. This is consistent with the observed age-redshift envelope, 
and it seems to indicate that we live in a universe with a classical vacuum energy density 
and that, on average, stars in the oldest galaxies formed about 0.7 Gyr after the Big Bang 
(e.g., at redshifts z ~ 7 — 10). To make the observed age-redshift relation consistent with 
a model with widely different w(z) behavior, we would have to infer that early-type galax- 
ies at different redshifts have very different formation epochs for their stellar populations. 
For example, the dotted line in Figure 4 indicates the age-redshift relation for an evolving 
population in a universe where w(z) = —2 for z > 1 and grows linearly from w = —2 at 
z = 1 to w = at z = 0. For this model to work, we would have to conclude that the oldest 
galaxies at z ~ 1.5 formed about 0.7 Gyr after the Big Bang, but that the oldest galaxies 
at z ~ formed ~ 3 Gyr after the Big Bang. Since the oldest galaxies in the LRG sample 
are as old as the oldest GCs, we consider this to be an unlikely explanation. Furthermore, 
not only would this scenario require that all local galaxies either formed their stars more 
recently than galaxies at high redshift or had a recent episodes of recent star formation, but 
would also require a remarkable fine-tuning for them all to obtain exactly the same age. We 
therefore conclude that this extreme model for w(z) is unlikely given the data. 

Unfortunately, the small number of galaxies in our samples at z > 0.2 does not allow us 
to compute dz/dt with enough accuracy to constrain w(z), and therefore this measurement 
will have to await better data. However, the z < 0.2 region is well populated by LRG 
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galaxies (Figure 4) and we can determine dt/dz at z ~ 0. This is illustrated in the bottom 
panel of Figure 4 by the clear shift in the upper envelope of the age histogram between the 
two redshift ranges < z < 0.04 and 0.08 < z < 0.12. In the next section we will investigate 
this phenomenon in more detail and use it to determine H . 

3.4. The value of H 

By concentrating on the age-redshift relation for SDSS LRG galaxies at z < 0.2, we 
now determine the Hubble constant Hq. This procedure relies upon determining the "edge" 
of the galaxy age distribution in different redshift intervals. 

If the ages of the galaxies were known with infinite accuracy, for each galaxy j at redshift 
Zi, one could associate an age Oj. The probability for the age of the oldest stellar population 
in galaxyj, P(t\a,i), would be given by a step function which jumps from to 1 at Oj. 
An age-redshift relation edge could then be obtained by dividing the galaxy sample into 
suitably-large redshift bins and multiplying the P{t\a,i) for all the galaxies in each bin. 

In practice, the age of each galaxy is measured with some error 5a,i. We thus assume 
P(t\a,i) = 1 if t > a, + 5 ai and lnP(aj) = — x 2 where x = (a* + 5 ai — t)/(y/25 ai ) otherwise. 
We have divided the z < 0.2 portion of the LRG sample into 51 redshift bins. For each 
bin, we obtain P(z,t) by multiplying the P(t\a,i) of the galaxies in that bin. We define 
the "edge" of the distribution t(z) to be where In P(z, t) drops by 0.5 from its maximum 
(i.e. AlnP = 0.5. We associate an error to this determination 5 t ( z ) given by t 2 (z) — t(z) 
where t 2 (z) corresponds to where In P(z,t) drops by 2 from its maximum (i.e. AlnP = 2, 
this approximately corresponds to the 68% confidence level). This procedure makes the 
determination of the "edge" less sensitive to the outliers (see also Raychaudhury et al. 
(1997)). Figure 6 presents the resultant t(z) relation. 

For z < z max , we fit t(z) with a straight line whose slope dt/dz is related to the Hubble 
constant at an effective redshift by H(z cS ) = — 1+ * (^~)~ 1 - This fit is performed by standard 
X 2 minimization. We also compute P> x a , the probability of obtaining equal or greater value 
of the reduced \ 2 if t ne t(z) points were truly lying on a straight line. Values of P> x a < 0.1 
means that a straight line is not a good fit to the points. 

For w = -1, we obtain H = H(z cS ) [Q m (l + ^ cff ) 3 + Q A V^ where Q m = 0.27, Q A = 0.7. 
Figure 7 shows how the H measurement depends on z max . For z max > 0.17 not only the 
value for Hq drifts but also a straight line is not a good fit to the points (P> x a suddenly 
drops below 0.1). For our H Q determination, z max = 0.17, P> x 2 = 0.32, z e s = 0.09, and the 
correction from H(z e s) to H is a 4% effect. We obtain H = 69 ± 12 km s _1 Mpc -1 . 
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Our H determination is in good agreement with the Hubble Key Project measurement 
(h = 0.72 ± 0.03 ± 0.07; Freedman et al. 2001), with the valued derived from the joint 
likelihood analysis of WMAP + 2dFGRS + Lyman-a forest power spectrum (h = 0.71^; 
Spergel et al. 2003), with gravitational lens time delay determinations {h = 0.591qq7 ±0.03; 
Treu & Koopmans 2002), and Sunyaev-Zeldovich measurements (h = O.GO^gi^is 13 ; Reese 
et al. 2002). We emphasize that our method is fully independent from any of the above 
measurements; it is reassuring that so many different methods yield such a consistent value 
for H Q . 



4. Conclusions 

We have obtained tight constraints on the effective value of w using solely the determina- 
tion of the ages of the oldest Galactic globular clusters and the observed position of the first 
acoustic peak of the CMB power spectrum. Our constraint on w is in good agreement with 
that of Spergel et al. (2003), but is subject to different possible systematics and is very weakly 
model dependent. Our constraint is slightly less stringent than that obtained by Spergel et 
al. (2003) from a joint likelihood analysis of WMAP with six external data sets ( WMAP + 
CBI + ACBAR + 2dFGRS + Lymana forest power spectrum + Type la supernovae + HST 
Key Project H determination), but is tighter than the CMB-only ( WMAP+CBI+ ACBAR) 
determination and comparable to the WMAP + ACBAR + CBI + HST constraint (Spergel 
et al. 2003). We speculate that the joint constraint from WMAP and all the external data 
sets considered in Spergel et al. (2003) could be further improved by adding stellar popula- 
tions age determination as an additional external data set. We have restricted our analysis 
to values of w > — 1 since theoretical models for w < — 1 seem difficult to construct (Carroll 
et al. 2003). The use of external data sets (see e.g. Kujat et al. (2002) and references therein 
for other approaches to measure w) could help to remove this restriction and constrain w 
for values < — 1. 

We have dated a large sample of galaxies and shown that a clear age-redshift envelope 
exists. By computing the relative ages of SDSS red galaxies at z < 0.2 we obtain an 
estimate of Hq independent of and in good agreement with determinations from the Hubble 
Key Project, gravitational lensing time delays, Sunyaev-Zeldovich, and WMAP. We also 
show that when surveys with a substantial number of high S/N spectra of old galaxies at 
z » become available, a determination of the possibly time- dependent equation of state 
parameter w(z) will be obtainable. 
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Fig. I. — The region of parameter space allowed if we fix 9a (which in turns fixes the position 
of the first CMB acoustic peak) and uo b = VL b h 2 = 0.02, weakly constraining 0.05 < VL m < 0.45 
and 0.5 < h< 0.9. 
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Fig. 2. — Confidence regions (68 and 90%) for the value of w are shown by the filled grey 
areas (marginalized over the position of the first acoustic peak in the CMB, The solid 
lines denote contours which show the constraints in w using the derived age of the oldest 
Galactic globular clusters only. The dash and dotted lines are the the position of the first 
acoustic peak and confidence region as measured by WMAP. Note that the position of the 
first peak greatly constrains the value of w. 
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Fig. 3. — The top 3 panels show examples of spectra from the SDSS LRG, Treu et al. and 
MS1054— 0321 samples. The dashed lines correspond to la error bars in the flux, the thin 
solid line is the measured spectrum and the thick solid line is the best fit to a single stellar 
population. The 3 bottom panels show confidence contours (68 and 95.4 %) in the recovered 
parameters. 
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Fig. 4. — Top panel: Age-redshift envelope obtained for the galaxies studied in this paper. 
Open circles correspond to galaxies from the SDSS LRG sample, triangles to the Treu et al. 
sample, diamonds to MS1054— 0321, and crosses are 53W091 and 53W069. A clear trend 
is present: galaxies age as the redshift decreases. The overall shape of this trend is in 
fair agreement with theoretical expectations for a ACDM (solid line) with Q m = 0.27 and 
H = 71 km s _1 Mpc -1 . The dotted line represent an alternative model disfavored by the 
data: w(z) = — 2 for z > 1 and then grows linearly from w = — 2 at z = 1 to w = at z = 0. 
Bottom panel: histogram along the age axis of top panel for the redshift range < z < 0.04 
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Fig. 5. — Simulated distribution for the age-redshift relation assuming that 10% of the mass 
was added to the early-type galaxy in the redshift interval < z < 1.5. We simulate two 
cases: a) there are only six early-type galaxies observed per age bin (top panel; similar 
to the higher-redshift samples), and b) there are 60 galaxies per age bin (bottom panel; 
corresponding to the SDSS sample). Note that the although the edge is less populated due 
to the fact of recent star formation, it is still clearly identifiable in both cases and is not 
biased. 
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Fig. 6. — The binned, low-redshift age-redshift relation, derived from the SDSS LRG sample. 
The solid line is a best-fit to the edge. 
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Fig. 7. — The value for H derived from Figure 4 as a function of redshift cut-off for com- 
puting the slope. Note that it is very weakly dependent on the value of this cut-off. The 
solid circle is our adopted value. The shaded are corresponds to the la confidence region 
from the Hubble Key Project. For z > 0.18 the probability of the fit being a straight line is 
P> x 2 < 0.1, while at z — 0.17 is 0.32. (see more details in text). 



